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Abstract. This paper presents a method for obtaining an analytic expression for the density 
function of observables in multifield models of inflation with sum-separable potentials. The 
most striking result is that the density function in general possesses a sharp peak and the 
location of this peak is only mildly sensitive to the distribution of initial conditions. A 
simple argument is given for why this result holds for a more general class of models than 
just those with sum-separable potentials and why for such models, it is possible to obtain 
robust predictions for observable quantities. As an example, the joint density function of the 
spectral index and running in double quadratic inflation is computed. For scales leaving the 
horizon 55 e-folds before the end of inflation, the density function peaks at n s = 0.967 and 
a = 0.0006 for the spectral index and running respectively. 
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1 Introduction 

The observable consequences of the simplest models of inflation with only one canonical 
field are quite well understood but ideas from fundamental physics seem to motivate models 
with more than one active field. 1 Models involving more than one light field, often referred 
to as "multifield models" allow for much richer inflationary behaviour and consequently 
making a prediction for observables in models of this kind is more complex. Although the 
phenomenology of these models is well studied, and the search for observational signatures 
such as local non-Gaussianity is an active area of research, at present it is not clear what 
even the simplest multifield models actually predict. 

One characteristic which makes computing observables more challenging is the fact that 
the primordial curvature perturbation evolves on superhorizon scales. In order to understand 
the possible signatures of multifield inflation it is essential to be able to precisely track the 
evolution after horizon crossing. Fortunately this is well studied and a number of techniques 
exist [2-28]. However, there is another important distinction between single field models of 
inflation and models with multiple active fields which is the sensitivity to initial conditions. 

Models involving only one field are essentially insensitive to initial conditions; some 
minimum number of e-folds of inflation is required to solve cosmological problems such as 
the flatness and homogeneity, but provided the total amount of inflation is a few e-folds more 
than this, the observational consequences are independent of how much inflation actually 
takes place. 2 The primordial curvature perturbation on a given scale is uniquely determined 

As described in Ref. [1], although it remains difficult in string theory to achieve a sufficiently light scalar 
to give rise to inflation, once this is managed, the mechanism which gives rise to one light scalar, tends to 
give rise to many. 

2 An important exception is the situation where there is more than one metastable vacuum. This situation 
certainly is sensitive to initial conditions. The discussion in this paper applies when only one minimum is 
relevant. 
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Figure 1. Sketch showing the horizon crossing surface in a single field inflationary potential and a 
two-field model. For a single field model, the horizon crossing surface is a single point, meaning that 
observables are approximately insensitive to inflationary dynamics prior to crossing this point. For 
the case of the two-field model, without knowledge of initial conditions one can not identify a single 
inflationary trajectory that corresponds to what we observe. One must therefore consider all possible 
inflationary trajectories. 

by the value of the inflaton at the moment that scale left the horizon. If we take the largest 
scales we observe to have left the horizon say 55 e-folds before the end of inflation, our 
observations are generally insensitive to what happened before then, be it another 20 e-folds 
of inflation or 200. 

When there is more than one active field the situation changes. As illustrated in Fig. 1, 
if the model involves Nf inflationary fields, instead of there being only one possible value 
of the inflaton at horizon crossing, now there is an infinite set of possible locations in field 
space forming an JVf — 1 dimentional hypersurface. 3 Without specifying initial conditions, 
it is not possible to say which inflationary trajectory corresponds to our observable universe 
and since different trajectories will in general give rise to different values for observables, the 
model is only as predictive as the volume in the space of observables permitted by the model. 
If we are to seriously confront multifield models of inflation with observation, then it is of 
paramount importance that this problem is overcome. 

The obvious question then is whether a description of initial conditions can be derived 
within the framework of the model. One can hope that an ultraviolet complete theory of 
inflation will provide information on the initial state but calculations along this line are 
clearly well beyond our current understanding of fundamental physics. A more promising 
approach is to consider how chaotic inflation populates the potential. Although this issue 
has received some attention in the past for the case of single field inflation (see for instance 

3 Discussed in more detail in section §2, this assumes only one trajectory passes through a given point in 
field space. 
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Refs. [29, 30] and references therein), a general description for the multifield case currently 
does not exist. 

Ultimately the most problematic aspect of this question is what is often referred to 
as the measure problem (see Ref. [31] for a recent review). Often described in the context 
of tunnelling between metastable vacua but nevertheless still of critical importance even to 
much simpler models [32] , the issue is that one must choose a measure when addressing what 
proportion of an infinite space-time corresponds to a particular choice of an infinite set of 
initial conditions. 

Despite the significant challenge of computing initial conditions, the aim of this paper 
is to show that at least for some models it is still possible to make robust predictions for 
observables. A general method is not given here. Instead, by considering the subclass of 
multifield models where there is no couplings between the fields, it is shown that the density 
function for observables can be computed analytically given an initial density function on 
field space. It turns out that under mild assumptions, the density function of observables is 
strongly peaked. This characteristic is largely determined by the geometry of the potential 
and can be comparatively insensitive to the initial density function on field space. This 
makes clear how the prediction for observables is (marginally) affected by the choice of 
measure, leading one to conclude that only certain choices will significantly change the result 
in comparison to the resolution of current and future observations such as Planck. 

The approach taken in this paper utilises the fact that for sum-separable potentials, 
provided isocurvature modes have decayed before the end of inflation, observables o may be 
expressed solely in terms of quantities evaluated at horizon crossing. As is often the case, 
in what follows this will be referred to as the horizon crossing approximation [33]. The use 
of this approximation means that all the information required to compute observables for 
all possible inflationary trajectories is contained in a single hypersurface £f parameterised 
by Nf — 1 variables {9\, . . , $/Vf _i } ■ The idea then is that by specifying the density function 
f{9\, ..,6>7v f _ 1 ) on the hypersurface, conservation of probability implies it is possible to com- 
pute the resulting density function for observables p(o\, o n ). For instance, consider the 
case where the hypersurface is 1-dimensional: Let 9 be a continuous random variable with 
probability density f(0). Defining o = o(9), provided the function is bijective, the probability 
density of o is given by p(o) where 

p(o)\do\ = f(8)\d6\. (1.1) 

The density function f(9) will depend on details of the theory as well as whatever the 
resolution of the measure problem may be, but even at this point Eq. (1.1) shows that a 
stationary point in o{6) will give rise to a spike in p(o) and so, with only mild assumptions 
as to the form of f(9), one can obtain a sharp prediction for the model. Furthermore, if the 
surface of evaluation is closed, then o(9) will be periodic and hence the existence of stationary 
points is guaranteed! 

With regard to relevant work existing in the literature, there seems to be little seeking 
to directly address the issue of predictions in multifield models. In Refs. [32, 34] a toy model 
of inflation in the landscape was investigated. By using a method equivalent to taking a flat 
distribution over initial conditions, distributions for observables were computed. However, 
the presence of multiple minima meant that the results were clearly sensitive to the choice 
of distribution of initial conditions yet it was not clear how to asses the impact of this choice 
on results. A similar approach was taken in Refs. [35-37] to study 6-field D-brane inflation 
but this model suffered from much the same problem. Here, drastically simpler models are 
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considered and the approach taken to studying these models is very different. The most 
important difference is that here the role of the initial density function is explicit. 

Despite not being directly about inflation, it is interesting to note the recent work of 
Sumitomo and Tye on the cosmological constant [38-40]. The mechanism by which they 
argue for the smallness of the cosmological constant is closely related to why a spike in the 
density function of observables is argued to be generic here. 

The rest of this paper is structured as follows. §2 briefly reviews the relevant expressions 
for observable quantities in sum-separable potentials and discusses the conditions under which 
the method sketched above is applicable. §3 is dedicated to a more thorough explanation of 
the use of Eq. (1.1). §4 is a detailed discussion of the prediction of quadratic inflation; as well 
as being an interesting model in its own right, the spherical symmetry of the model makes 
calculating predictions exceptionally straightforward and hence will serve as a transparent 
demonstration of the method. §5 briefly discusses why the results should hold for a more 
general class of models than just the class of sum-separable models and §6 concludes this 
paper. 

2 Sum-Separable Potentials and the Horizon Crossing Approximation 

This section briefly introduces the expressions for observables in canonical models of inflation 
with sum-separable potentials and then discusses the conditions necessary for the proposed 
method to be applicable. 

2.1 Expressions for Observables 

The curvature perturbation may be related to perturbations in the fields by realising that 
on large scales £ is equivalent to the perturbation of the number of e-foldings from an initial 
flat hypersurface at t = t*, to a final uniform-density hypersurface at t = t c [2-4, 8] 

C(t c , x) ~ 5N(t c , U,x) = N{t c , U,x)- N(t , U). (2.1) 

where N(t c , t*) = f£ Hdt. Expanding 8N in terms of the initial field perturbations to second 
order, one obtains 

C(t c , x) = SN{t c , U,x) = N, a 5<f>* a + ^N, a p 54>*J4>% (2.2) 

where repeated indices should be summed over, and N, a , N, a p represent first and second 
derivatives of the number of e-folds with respect to the fields </>* . 

Observables of interest are related to the correlation functions of £ and hence may be 
expressed in terms of field correlation functions multiplied by derivatives of N . In general 
computing the N derivatives require numerical techniques but a useful exception is the case 
of sum-separable potentials 

Nt 

W{cj> l ,..,4> Ni )=Y J V i {cj> i ). (2.3) 

i 

Models of this kind have the appealing characteristic that observables such as the power 
spectrum, spectral index, running, non-Gaussianity parameter /nl etc. may be computed 
analytically. It should be noted however that there is no particularly good reason to believe 
that these models are well motivated from an effective field theory perspective; one in general 
should expect couplings between the fields. Although a more general approach is not given 
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here, in later sections it will be argued that the main result of this paper applies to a broader 
class of models than just those with sum-separable potentials. 

This paper focusses on quantities relating to the two-point statistics of £, particularly 
the spectral index n s and running a, however all techniques used throughout this work are 
just as easily applied to higher order statistics. The reader is referred to Refs. [36, 41, 42] 
for derivations of the expressions. For this paper it will suffice to quote the results. Setting 
Mpl = 1, the power spectrum is given by 

rxr n f 2 



e 

a=l c 



where 



— ^ (2.5) 



and e a is the first slow-roll parameter 



1 (V^ 2 



(2.6) 



2 \W, 

such that e = ^ e a . Also, rather importantly for what follows, the term 

^ = -EW-e c M (2-7) 

contains all the information about the constant density surface at the time of evaluation and 
all other terms are evaluated at horizon-crossing. 

By differentiating with respect to Ink, the spectral index is found to be 

*.-l = -2e*-4 V Np g^A (2-8) 
where r) a is the second slow-roll parameter 

V" 

Va - (2-9) 
Differentiating with respect to In k once again, the running is found to be [36] 

a = -8e* 2 +4Ve^-16^ 8 ^ ^+4e*^^f 2 \ 

a=l O a # J ^ ^ ^ a l—jol e * 



(2.10) 

where £ Q is the third slow-roll parameter 

V'V" 

(2.11) 

All of the above expressions can be shown to reduce to the standard single-field formula by 
setting u a = 1. 
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2.2 The horizon crossing approximation and the adiabatic limit 

As already discussed, the method proposed here for mapping a density function in field space 
to a density function for observables makes use of the horizon crossing approximation. The 
above expressions for the power spectrum (2.4), tilt (2.8) and running (2.10), are all composed 
of slow-roll terms evaluated at horizon crossing "*" and one other term u a , of which the only 
contribution not evaluated on "*" is Z a . It follows that in order for this approach to be 
appropriate two requirements must be satisfied: 

1. The slow-roll approximations e <C 1, )] < 1, ( <C 1 must be valid. 

2. Observables must stop evolving within the realm of validity of the model. 

The first requirement is necessary for a number of reasons. The method proposed here seeks 
to map the density function on an JVf - 1 hypersurface in field space to a density function for 
observables. The very fact that the objective is phrased this way already implicitly assumes 
that (up to a time shift) there is a unique trajectory passing through each point in field 
space. A simple way to ensure this is the case is to stipulate that the slow-roll conditions 
are satisfied. If one does not wish to place any constraints on momenta then the situation 
is drastically more complicated. In principle an infinite number of both inflationary and 
non-inflationary trajectories may pass through any given point in field space. To handle such 
a situation is well beyond the scope of this approach. 

Another reason for the first requirement is that the expressions for observables given 
above are obtained by assuming slow-roll. In practice it can well be the case that these 
expressions are valid even when the slow-roll approximations start to break down but this 
range of validity is clearly model-dependant. 

The second requirement refers to the fact that for the model to be predictive, the primor- 
dial curvature perturbation must become conserved before the end of inflation. Otherwise 
a description of reheating is required. In order for the curvature perturbation to become 
conserved, isocurvature modes must have exponentially decayed by the end of inflation such 
that the model approaches the adiabatic limit [32, 36, 43-45]. Although it is difficult to say 
with certainty when isocurvature modes are sufficiently decayed, an intuitive method to test 
that this requirement has been satisfied is to track the evolution of the width of the bundle 
of perturbed inflationary trajectories in field space. The reader is referred to Ref. [45] for 
detailed discussion but an expression for the bundle width can be found to be 



where the expression given here differs from that given in Ref. [45] slightly since this has been 
written for the case of sum-separable potentials and made use of dN = — J2 a Va/V^d^a. For 
an adiabatic limit to be reached, a necessary (but not sufficient if N{ > 2) condition is that 
O — > 0. By inspection of Eq. (2.12), a period of focusing requires 



though how much focussing is required is a model dependant statement. 

Satisfaction of the second requirement guarantees Z a will be exponentially close to a 
constant at the end of inflation but does not guarantee that it is negligible. Satisfaction of 




(2.12) 




(2.13) 



a 
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both the first and the second requirement does however guarantee that to a good approxima- 
tion one can express observables solely in terms of quantities evaluated at horizon crossing. 
This is what is referred to as the horizon crossing approximation. 4 In practice, it is often 
the case that Z a — > in the approach to the adiabatic limit but even in the situations where 
this is not the case, provided the first requirement is satisfied, it will be possible to map a 
given location on the final surface to a unique location on the horizon crossing surface. For 
the example of quadratic inflation given in this paper Z a — > at the end of inflation. 



3 Computing the density function of observables 

The objective is to obtain an expression for p(o\, .., o n ) where Oj are the observables of interest 
such as the spectral index, running etc. If the model consists of Nf fields, then the horizon 
crossing contour is an (iVf — l)-dimensional hypersurface in field space and hence, under 
the horizon crossing approximation, each observable may be expressed in terms of Nf — 1 
parameters 6 a i. Latin indices label the observables, greek indices label coordinates on field 
space and primed greek indices label coordinates on the horizon crossing hypersurface in 
field space. The choice of parameterisation is arbitrary. Describing the horizon crossing 
hypersurface parametrically such that <j) a = cp a (9i, .., 0n { _ 1 ), the volume element on the 
horizon crossing hypersurface can be expressed in terms of the parameters 9 a i in the usual 
way. Constructing basis vectors 

e *' ~ de a ,> (3 - 1} 

assuming flat field space, the induced metric is then 



so remembering that the volume of the corresponding parallelepiped is given by \J det g Sf , 
the volume element d9 in Eq. (1.1) can now be written more explicitly as 



d6 -»■ Vdet g^d9 x ■ ■ ■ dB Nt -i, (3.3) 

If n < Nf then a similar picture holds for the left hand side of Eq. (1.1). Oj(#i, .., #jv f _i) 
are the parametric equations of a surface in the space of observables and so one can construct 
the volume element in the same way. The induced metric in "observables space" is 

>A doi doi 



and so Eq. (1.1) becomes 



/ det (? Sf 

p(oi,..,o n ) = f[(j)x (#!,.. ,0jv f -i),..,<?!>jv f (#1, ..,0jv f -i)]y det gSo ■ (3-5) 

Although it will always be possible to compute this expression when the horizon crossing 
approximation is valid and n < Nf, clearly some sacrifices are made to achieve this generality. 

4 N.B this is not the same as evaluating the expression for a given observable at horizon crossing! 
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Ideally what one would like to do is perform the inverse mapping and express 8 a '(oi, o n ) 
such that p(oi, ..,o n ) can be written explicitly in terms of the observables. Whether this is 
practical is a model dependant statement and the difficulty in doing so will also depend on 
the number of observables under consideration. For instance if n > Nf, then one either needs 
to compute the cumulative distribution function or construct enough dummy variables such 
that n = Nf. To give an explicit example, consider the case of 1 observable n = 1 and a 
2-dimensional horizon crossing surface Nf — 1 = 2. One would essentially need to compute 

p(o) = Y jj fVd^d9fde 2 , (3.6) 

which in general one probably would prefer to avoid in favour of numerical techniques. 

However Eq. (3.5) does have some advantages. As already mentioned, it is simple to 
compute and still allows one to locate and compare peaks in the density function. In some 
cases, it is sufficient to construct the full density function and even make a sensible estimate 
of confidence limits. The main advantage as far as this paper is concerned is that Eq. (3.5) 
makes manifest the relationship between the geometry of the potential and the initial density 
function. One potentially valuable application of this is in understanding the sensitivity of 
the model to the choice of measure. This will be discussed in the context of the example in 
the following section. 



4 Quadratic Inflation 

In this section an explicit example is given of the method outlined in §3. Double quadratic 
inflation is a good example for a number of reasons. Explicitly shown in Ref. [45] for an 
arbitrary number of fields approaching a quadratic minimum, all isocurvature modes are 
exponentially decaying before the end of inflation and so the horizon crossing approximation 
is valid. The spherical symmetry of the horizon crossing surface also makes it a particularly 
simple example. 

The potential for double quadratic inflation is [46] 

Vi = V 2 = \m 2 2 4>l (4.1) 

such that W = V\ + V 2 . An expression for the horizon crossing contour may be obtained by 
writing the number of e-folds as 



and so 




A well known [42, 46] and rather helpful choice of parameterisation is to move into polar 
coordinates. By neglecting the contribution from the "c" surface one can write 

01 = 2VNcos9 (f> 2 = 2\/iVsin6», (4.4) 
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Figure 2. Example plots of n s {9) and a{9) for scales leaving the horizon 55 e-folds before the end of 
inflation in double quadratic inflation with masses 777.2/mi = 9. 



where here and what follows, unless explicitly written otherwise, fields only lie on the surface 
Sf and so the "*" label has been dropped. For quadratic inflation, the observables (2.4), 
(2.8) and (2.10) become 



r c (e,N) 

n s (6,N) 

a(6,N) 



1 



El 

4tt 2 
-2e 

-8e 2 



1 

N' 
2 



+ 4(ei?7i + e 2 m) ■ 



(4.5) 
(4.6) 
(4.7) 



The amplitude of the power spectrum is in a sense less interesting than the other observables 
since it may always be adjusted by a pre- factor on the potential which does not affect the 
inflationary dynamics. Other observables one might want to consider are the non-Gaussianity 
parameter /nl and the tensor-to-scalar ratio r. For the case of double quadratic inflation, it 
was shown in Ref. [42] that -§/nl = 1/N. Since £f is defined as being at a fixed number 
of e-folds before the end of inflation, /nl is single valued over T,f which is why it has not 
been part of the discussion until now (the same is true for the tensor-to-scalar ratio). Hence, 
for this example, the objective is to calculate p(n s ,a). Fig. 2 shows plots of the spectral 
index and running evaluated 55 e-folds before the end of inflation. As mentioned previously, 
an important characteristic to bear in mind for what follows is that these observables are 
periodic in 9. This is simply a consequence of the fact that represents coordinates on a 
closed surface. It is also worth noting that the stationary points at = ni/2, where i is an 
integer, occur at the same values for n s and a. 

4.1 Regarding the field space density function fs f 

It is necessary to specify the density function /s f . As already mentioned, this is dependant 
on the details of the model as well as the choice of measure. For the purposes of this paper 
the distribution is chosen to be flat 

1 



(4.8) 



Following Ref. [47] , this choice is a statement of ignorance. We simply adopt the distribution 
requiring the least additional assumptions. No further justification will be given at this stage, 
however in §5.1 the implications of this choice will be discussed. 
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Figure 3. Example plots for double quadratic potential with masses mijmx = 9 (red), TO2/W1 = 7 
(orange), m-xjrax = 5 (gold). The plot on the left shows the set of possible values of {n s ,a}. The 
right hand plot additionally includes p(n s ,a). The density function for all examples shown peak at 
n s — 1 = — and a = — and are divergent and so the full plot range is not shown. The grey 
lines are the projections onto the p{n s , a) — plane. 



The physics in each quadrant of field space is the same, so without loss of generality 
consider just the first quadrant where both fields are positive. The density function over the 
contour is then 

h - = WV <4 - 9) 

4.2 The joint density function p(n s ,a) 

As mentioned, it is convenient to parameterise in terms of polar coordinates {N, 9}. Pythagoras 
and Eq. (4.4) give cEf = 2^/NdQ and so the right hand side of Eq. (1.1) becomes 

fedZt = -d9. (4.10) 

Since for this model the observables of interest are o = {n s , a}, using Eq. (3.1) and Eq. (3.2) 
the element do may be rewritten as 




do = \H^) +{^) 2 de (4.11) 
and hence combing Eq. (1.1) with Eq. (4.10) and Eq. (4.11), 

p(n s ,a)= 2 ^ 1 (4.12) 

*7(*) +($) 

The right hand plot of Fig. 3 shows p(n s ,a) for the case of double quadratic inflation 
with masses 7712/mi = 9, mil mi = 7, mijm\ = 5 shown in red, orange and gold respectively. 
The left plot shows the contours of possible values in the n s — a plane. There is a strong peak 
in the density function which corresponds 9 = iri/2, where % is a member of the integers. 
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Comparing with Fig. 2, this is exactly what should be expected. Fig. 2 shows that for a large 
range of 9, the tilt and running are slowly varying. Since a large proportion of the contour 
maps to a relatively small proportion of the possible range of values of the tilt and running, 
this will give rise to a very sharp spike in the density function. 

It is interesting to note that the peak occurs at the same point for each mass ratio, i.e 
the prediction of the model is not as sensitive to the mass hierarchy as the space of possible 
values for observables might suggest. 

5 Discussion 

In computing the density function for observables in the previous section, two important 
assumptions were made. One was that the density function on field space was taken to be 
a flat distribution over the horizon crossing contour. The second was that the potential was 
sum-separable. In this section the significance of these assumptions is assessed. 

5.1 The choice of a flat density function at horizon crossing 

Following the procedure in Ref. [46], an expression for the evolution of 9 from some initial 
surface £j to the horizon crossing surface S* can be obtained by substituting Eq. (4.4) into 
the slow-roll equations of motion for the fields and integrating 



As shown by the example plots in Fig. 4, there is a dynamical attractor causing the 
field trajectories to converge on 9 = or 9 = tt. This means that a flat distribution chosen 
on some N > N* contour will evolve to give rise to a density function at horizon crossing 
with peaks at 9 = and 9 = n. As shown in Fig. 4, the result of this is to sharpen the peak 
in the density function of observables. In this sense choosing a flat distribution at horizon 
crossing can be considered to result in a lower bound on the strength of the peak sourced by 
dynamical effects. 

As mentioned previously, one advantage of computing the density function as given by 
Eq. (3.5), is that it makes clear the role of the density function on field space and so in turn, 
the choice of measure. Since field space dynamics act to strengthen the peak in p(o), in order 
for the prediction of the model to be changed significantly, one would require a choice of 
measure that acts to counter this dynamical effect. 

5.2 The use of sum-separable potentials 

Central to the approach taken in this paper is the use of sum-separable potentials. It is 
therefore important to understand to what extent the results are special to this class of 
models. The reason the analysis performed here was restricted to this class was that it 
enabled the use of the horizon crossing approximation and hence provided a simple method 
of mapping a density function on field space to a density function for observables. The most 
striking result is that the functional form of Eq. (3.5) seems to imply that the density function 
of observables will be strongly peaked and this is indeed the case for the example of double 
quadratic inflation that was explored. If a strong peak turns out to be a generic feature 
of multifield models this is advantageous on two fronts. Most importantly it means that a 
given multifield model may be significantly more predictive than considering just the range 




(5.1) 
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Figure 4. Example plots showing the effect on p(n s — 1, a) if a flat distribution is taken at 200 e-folds 
before the end of inflation instead of at horizon crossing, for the case of double quadratic inflation 
with mass ratio m^/mi = 9. The top left plot shows the evolution of example trajectories from 200 
e-folds before the end of inflation up until horizon crossing. The plot shows a significant proportion 
of the trajectories converge at 9 = and = %. The top right plot shows the relation of 9i, the 
parameterisation of the initial contour, to 0„, the parameterisation of the horizon crossing contour. 
The plot shows once again that trajectories over a significant range of Oi converge onto 0* = and 
9* = 7T. The bottom left hand plot shows the resulting density function over one quadrant of the 
horizon crossing surface with a strong peak at 9 = 0. Finally, the bottom right plot shows the result 
on p(n s — 1, a) in gold compared to p(n s — 1, a) for a flat distribution at horizon crossing shown in 
red. The red plot is the same as that shown in the right hand plot of Fig. 3 but projected onto the 
p(n s — 1, a) — n s — 1 plane. The effect of the evolution is to sharpen the peak in p(n s — 1, a) — n s — 1. 



of possible values of observables allowed within the model and hence one can be much more 
optimistic about the prospect of constraining such models. Secondly, if sharp peaks in the 
density function exist, then knowing the location of these peaks can be very advantageous 
to numerical approaches to computing the density function. 

Although the use of the horizon crossing approximation provides a simple map from field 
space to observables, it is argued here that it is not the cause of the peak in p(o). It may not 
be very efficient but clearly there already exists a method of computing the density function 
of observables for a much broader class inflationary models than the method given here. For 
example, consider the class of all inflationary models where the slow-roll approximation is 
valid at horizon crossing and also reach an adiabatic limit before the end of inflation. 5 For 



5 The method proposed here is straightforward to adapt to other models of inflation but some details may 



o(6) 




'/(*) 

Figure 5. Sketch of why a stationary point in the functional form of 0(8) will give rise to a divergence 
in the density function p(o). The spacing of the red lines represents the density function f(ff) and 
how it is distorted under a change of variable to give p(o). 

this class one could perform the following procedure: 

1. Perform a Monte Carlo search for the horizon crossing surface £*. 

2. Draw initial field space positions from a given density function /j^. 

3. For each set of initial field space positions, starting with momenta given by the slow-roll 
equations, evolve each trajectory until the end of inflation. 

4. For each inflationary trajectory, using one of the methods given in Refs. [2-28], compute 
relevant observables. 

In doing this, one would numerically obtain a map from the iVf — 1 dimensional surface 
on which inflation ends to observables which may, as before, be parameterised by Nf — 1 
variables {0\, .., 9n {1 }. Once this map is obtained, Eq. (3.5) is once again applicable. 

If the surface in field space is closed there are only two options for the mapping. Either 
the observable is independent of {9±, .., #jv f _i}j or it is periodic in {9±, .., #jv f _i}- If the 
observable is independent of {9%, .., 0jv f _i} then (as was the case here for /nl and r) that 
observable is single valued. If the observable is periodic in {9\, .., ^Ar f _ x }, then stationary 
points are guaranteed. 

If the surface in field space is not closed then it is not possible to say anything so 
concrete. Nevertheless, the conditions for a peak in the density function to occur are clear. 
As illustrated in Fig. 5, a peak in p(o) occurs when do/d9 is small. If the surface in field space 
is not closed, then this requirement can not be guaranteed but it is certainly still permitted. 

vary and so in pursuit of being a little more concrete, discussion has been restricted to this class. 



- 13 - 



6 Conclusion 



If multifield models of inflation are to be confronted with data, it is essential that we un- 
derstand what the prediction of a given model actually is. As has been discussed, there is 
essentially two parts to this problem: 

1. Computing the density function on field space. 

2. Mapping the density function on field space to the density function of observables. 

Progress with the first part is stymied by the measure problem. This issue has received 
considerable attention recently and a number of authors have proposed various possible so- 
lutions. At present however, the problem remains open. The second part has received com- 
paratively little attention but it is hoped that this paper represents steps towards rectifying 
this. 

As an attempt at making progress with the second issue, here a solution for the case 
of sum-separable potentials was proposed. Although it is useful to have analytic examples, 
clearly this is not sufficient and a more general approach is urgently needed if we are to 
compare multifield models with the soon to arrive data from Planck. As discussed, in prin- 
ciple a numerical approach is already available. However it is inefficient, meaning that a 
thorough analysis of models with a large number of active fields, or models with complicated 
Lagrangians may still be computationally unaccessible. 

The main point this work is trying to make is that a strongly peaked density function 
for observables is expected to be generic for a broad class of models. This means that some 
models may be relatively insensitive to the density function on field space and hence for such 
models it should be possible to make robust predictions with the knowledge we already have 
available to us. 
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